Differential gene expression in breast cancer cell lines and stroma-tumor differences in microdissected breast cancer biopsies revealed by display array analysis

To examine gene-expression patterning in late-stage breast cancer biopsies, we used a microdissection technique to separate tumor from the surrounding breast tissue or stroma. A DD-PCR protocol was then used to amplify expressed products, which were resolved using PAGE and used as probe to hybridize with representative human arrays and cDNA libraries. The probe derived from the tumor–stroma comparison was hybridized with a gene array and an arrayed cDNA library derived from a GCT of bone; 21 known genes or expressed sequence tags were detected, of which 17 showed differential expression. These included factors associated with epithelial to mesenchymal transition (vimentin), the cargo selection protein (TIP47) and the signal transducer and activator of transcription (STAT3). Northern blot analysis was used to confirm those genes also expressed by representative breast cancer cell lines. Notably, 6 genes of unknown function were restricted to tumor while the majority of stroma-associated genes were known. When applied to transformed breast cancer cell lines (MDA-MB-435 and T47D) that are known to have different metastatic potential, DD array analysis revealed a further 20 genes; 17 of these genes showed differential expression. Use of microdissection and the DD-PCR array protocol allowed us to identify factors whose localized expression within the breast may play a role in abnormal breast development or breast carcinogenesis.

Design, construction, and analysis of cell line arrays and tissue microarrays for gene expression analysis

Cell line array (CMA) and tissue microarray (TMA) technologies are high-throughput methods for analysing both the abundance and distribution of gene expression in a panel of cell lines or multiple tissue specimens in an efficient and cost-effective manner. The process is based on Kononen's method of extracting a cylindrical core of paraffin-embedded donor tissue and inserting it into a recipient paraffin block. Donor tissue from surgically resected paraffin-embedded tissue blocks, frozen needle biopsies or cell line pellets can all be arrayed in the recipient block. The representative area of interest is identified and circled on a haematoxylin and eosin (H&E)-stained section of the donor block. Using a predesigned map showing a precise spacing pattern, a high density array of up to 1,000 cores of cell pellets and/or donor tissue can be embedded into the recipient block using a tissue arrayer from Beecher Instruments. Depending on the depth of the cell line/tissue removed from the donor block 100-300 consecutive sections can be cut from each CMA/TMA block. Sections can be stained for in situ detection of protein, DNA or RNA targets using immunohistochemistry (IHC), fluorescent in situ hybridisation (FISH) or mRNA in situ hybridisation (RNA-ISH), respectively. This chapter provides detailed methods for CMA/TMA design, construction and analysis with in-depth notes on all technical aspects including tips to deal with common pitfalls the user may encounter. © Springer Science+Business Media, LLC 2011.

Functional genes and gene array analysis as tools for monitoring hydrocarbon biodegradation

Bioremediation, which is the exploitation of the intrinsic ability of environmental microbes to degrade and remove harmful compounds from nature, is considered to be an environmentally sustainable and cost-effective means for environmental clean-up. However, a comprehensive understanding of the biodegradation potential of microbial communities and their response to decontamination measures is required for the effective management of bioremediation processes. In this thesis, the potential to use hydrocarbon-degradative genes as indicators of aerobic hydrocarbon biodegradation was investigated. Small-scale functional gene macro- and microarrays targeting aliphatic, monoaromatic and low molecular weight polyaromatic hydrocarbon biodegradation were developed in order to simultaneously monitor the biodegradation of mixtures of hydrocarbons. The validity of the array analysis in monitoring hydrocarbon biodegradation was evaluated in microcosm studies and field-scale bioremediation processes by comparing the hybridization signal intensities to hydrocarbon mineralization, real-time polymerase chain reaction (PCR), dot blot hybridization and both chemical and microbiological monitoring data. The results obtained by real-time PCR, dot blot hybridization and gene array analysis were in good agreement with hydrocarbon biodegradation in laboratory-scale microcosms. Mineralization of several hydrocarbons could be monitored simultaneously using gene array analysis. In the field-scale bioremediation processes, the detection and enumeration of hydrocarbon-degradative genes provided important additional information for process optimization and design. In creosote-contaminated groundwater, gene array analysis demonstrated that the aerobic biodegradation potential that was present at the site, but restrained under the oxygen-limited conditions, could be successfully stimulated with aeration and nutrient infiltration. During ex situ bioremediation of diesel oil- and lubrication oil-contaminated soil, the functional gene array analysis revealed inefficient hydrocarbon biodegradation, caused by poor aeration during composting. The functional gene array specifically detected upper and lower biodegradation pathways required for complete mineralization of hydrocarbons. Bacteria representing 1 % of the microbial community could be detected without prior PCR amplification. Molecular biological monitoring methods based on functional genes provide powerful tools for the development of more efficient remediation processes. The parallel detection of several functional genes using functional gene array analysis is an especially promising tool for monitoring the biodegradation of mixtures of hydrocarbons.

Lesion Detection in Breast Ultrasound Images Using Tissue Transition Analysis

Breast cancer is one of the leading cause of cancer related deaths in women and early detection is crucial for reducing mortality rates. In this paper, we present a novel and fully automated approach based on tissue transition analysis for lesion detection in breast ultrasound images. Every candidate pixel is classified as belonging to the lesion boundary, lesion interior or normal tissue based on its descriptor value. The tissue transitions are modeled using a Markov chain to estimate the likelihood of a candidate lesion region. Experimental evaluation on a clinical dataset of 135 images show that the proposed approach can achieve high sensitivity (95 %) with modest (3) false positives per image. The approach achieves very similar results (94 % for 3 false positives) on a completely different clinical dataset of 159 images without retraining, highlighting the robustness of the approach.

Teleseismic array analysis of upper mantle compressional velocity structure

Large quantities of teleseismic short-period seismograms recorded at SCARLET provide travel time, apparent velocity and waveform data for study of upper mantle compressional velocity structure. Relative array analysis of arrival times from distant (30° < Δ < 95°) earthquakes at all azimuths constrains lateral velocity variations beneath southern California. We compare dT/dΔ back azimuth and averaged arrival time estimates from the entire network for 154 events to the same parameters derived from small subsets of SCARLET. Patterns of mislocation vectors for over 100 overlapping subarrays delimit the spatial extent of an east-west striking, high-velocity anomaly beneath the Transverse Ranges. Thin lens analysis of the averaged arrival time differences, called 'net delay' data, requires the mean depth of the corresponding lens to be more than 100 km. Our results are consistent with the PKP-delay times of Hadley and Kanamori (1977), who first proposed the high-velocity feature, but we place the anomalous material at substantially greater depths than their 40-100 km estimate.

Detailed analysis of travel time, ray parameter and waveform data from 29 events occurring in the distance range 9° to 40° reveals the upper mantle structure beneath an oceanic ridge to depths of over 900 km. More than 1400 digital seismograms from earthquakes in Mexico and Central America yield 1753 travel times and 58 dT/dΔ measurements as well as high-quality, stable waveforms for investigation of the deep structure of the Gulf of California. The result of a travel time inversion with the tau method (Bessonova et al., 1976) is adjusted to fit the p(Δ) data, then further refined by incorporation of relative amplitude information through synthetic seismogram modeling. The application of a modified wave field continuation method (Clayton and McMechan, 1981) to the data with the final model confirms that GCA is consistent with the entire data set and also provides an estimate of the data resolution in velocity-depth space. We discover that the upper mantle under this spreading center has anomalously slow velocities to depths of 350 km, and place new constraints on the shape of the 660 km discontinuity.

Seismograms from 22 earthquakes along the northeast Pacific rim recorded in southern California form the data set for a comparative investigation of the upper mantle beneath the Cascade Ranges-Juan de Fuca region, an ocean-continent transit ion. These data consist of 853 seismograms (6° < Δ < 42°) which produce 1068 travel times and 40 ray parameter estimates. We use the spreading center model initially in synthetic seismogram modeling, and perturb GCA until the Cascade Ranges data are matched. Wave field continuation of both data sets with a common reference model confirms that real differences exist between the two suites of seismograms, implying lateral variation in the upper mantle. The ocean-continent transition model, CJF, features velocities from 200 and 350 km that are intermediate between GCA and T7 (Burdick and Helmberger, 1978), a model for the inland western United States. Models of continental shield regions (e.g., King and Calcagnile, 1976) have higher velocities in this depth range, but all four model types are similar below 400 km. This variation in rate of velocity increase with tectonic regime suggests an inverse relationship between velocity gradient and lithospheric age above 400 km depth.

Proteomic analysis of colorectal cancer metastasis:stathmin-1 revealed as a player in cancer cell migration and prognostic marker

Metastasis accounts largely for the high mortality rate of colorectal cancer (CRC) patients. In this study, we performed comparative proteome analysis of primary CRC cell lines HCT-116 and its metastatic derivative E1 using 2-D DIGE. We identified 74 differentially expressed proteins, many of which function in transcription, translation, angiogenesis signal transduction, or cytoskeletal remodeling pathways, which are indispensable cellular processes involved in the metastatic cascade. Among these proteins, stathmin-1 (STMN1) was found to be highly up-regulated in E1 as compared to HCT-116 and was thus selected for further functional studies. Our results showed that perturbations in STMN1 levels resulted in significant changes in cell migration, invasion, adhesion, and colony formation. We further showed that the differential expression of STMN1 correlated with the cells' metastatic potential in other paradigms of CRC models. Using immunohistochemistry, we also showed that STMN1 was highly expressed in colorectal primary tumors and metastatic tissues as compared to the adjacent normal colorectal tissues. Furthermore, we also showed via tissue microarray analyses of 324 CRC tissues and Kaplan-Meier survival plot that CRC patients with higher expression of STMN1 have poorer prognosis.

Whole genome methylation array analysis reveals new aspects in Balkan endemic nephropathy etiology

Background Balkan endemic nephropathy (BEN) represents a chronic progressive interstitial nephritis in striking correlation with uroepithelial tumours of the upper urinary tract. The disease has endemic distribution in the Danube river regions in several Balkan countries. DNA methylation is a primary epigenetic modification that is involved in major processes such as cancer, genomic imprinting, gene silencing, etc. The significance of CpG island methylation status in normal development, cell differentiation and gene expression is widely recognized, although still stays poorly understood. Methods We performed whole genome DNA methylation array analysis on DNA pool samples from peripheral blood from 159 affected individuals and 170 healthy individuals. This technique allowed us to determine the methylation status of 27 627 CpG islands throughout the whole genome in healthy controls and BEN patients. Thus we obtained the methylation profile of BEN patients from Bulgarian and Serbian endemic regions. Results Using specifically developed software we compared the methylation profiles of BEN patients and corresponding controls and revealed the differently methylated regions. We then compared the DMRs between all patient-control pairs to determine common changes in the epigenetic profiles. SEC61G, IL17RA, HDAC11 proved to be differently methylated throughout all patient-control pairs. The CpG islands of all 3 genes were hypomethylated compared to controls. This suggests that dysregulation of these genes involved in immunological response could be a common mechanism in BEN pathogenesis in both endemic regions and in both genders. Conclusion Our data propose a new hypothesis that immunologic dysregulation has a place in BEN etiopathogenesis. Keywords: Epigenetics; Whole genome array analysis; Balkan endemic nephropathy

High-sensitivity array analysis of gene expression for the early detection of disseminated breast tumor cells in peripheral blood

Early detection is an effective means of reducing cancer mortality. Here, we describe a highly sensitive high-throughput screen that can identify panels of markers for the early detection of solid tumor cells disseminated in peripheral blood. The method is a two-step combination of differential display and high-sensitivity cDNA arrays. In a primary screen, differential display identified 170 candidate marker genes differentially expressed between breast tumor cells and normal breast epithelial cells. In a secondary screen, high-sensitivity arrays assessed expression levels of these genes in 48 blood samples, 22 from healthy volunteers and 26 from breast cancer patients. Cluster analysis identified a group of 12 genes that were elevated in the blood of cancer patients. Permutation analysis of individual genes defined five core genes (P ≤ 0.05, permax test). As a group, the 12 genes generally distinguished accurately between healthy volunteers and patients with breast cancer. Mean expression levels of the 12 genes were elevated in 77% (10 of 13) untreated invasive cancer patients, whereas cluster analysis correctly classified volunteers and patients (P = 0.0022, Fisher's exact test). Quantitative real-time PCR confirmed array results and indicated that the sensitivity of the assay (1:2 × 108 transcripts) was sufficient to detect disseminated solid tumor cells in blood. Expression-based blood assays developed with the screening approach described here have the potential to detect and classify solid tumor cells originating from virtually any primary site in the body.

An orchestrated gene expression component of neuronal programmed cell death revealed by cDNA array analysis

Programmed cell death (PCD) during neuronal development and disease has been shown to require de novo RNA synthesis. However, the time course and regulation of target genes is poorly understood. By using a brain-biased array of over 7,500 cDNAs, we profiled this gene expression component of PCD in cerebellar granule neurons challenged separately by potassium withdrawal, combined potassium and serum withdrawal, and kainic acid administration. We found that hundreds of genes were significantly regulated in discreet waves including known genes whose protein products are involved in PCD. A restricted set of genes was regulated by all models, providing evidence that signals inducing PCD can regulate large assemblages of genes (of which a restricted subset may be shared in multiple pathways).

RUNX3 downregulation in human lung adenocarcinoma is independent of p53, EGFR or KRAS status

RUNX3 aberrations play a pivotal role in the oncogenesis of breast, gastric, colon, skin and lung tissues. The aim of this study was to characterize further the expression of RUNX3 in lung cancers. To achieve this, a lung cancer tissue microarray (TMA), frozen lung cancer tissues and lung cell lines were examined for RUNX3 expression by immunohistochemistry, while the TMA was also examined for EGFR and p53 expression. RUNX3 promoter methylation status, and EGFR and KRAS mutation status were also investigated. Inactivation of RUNX3 was observed in 70% of the adenocarcinoma samples, and this was associated with promoter hypermethylation but not biased to EGFR/KRAS mutations. Our results suggest a central role of RUNX3 downregulation in pulmonary adenocarcinoma, which may not be dependent of other established cancer-causing pathways and may have important diagnostic and screening implications.

The prognostic significance of the aberrant extremes of p53 immunophenotypes in breast cancer

Aims: The utility of p53 as a prognostic assay has been elusive. The aims of this study were to describe a novel, reproducible scoring system and assess the relationship between differential p53 immunohistochemistry (IHC) expression patterns, TP53 mutation status and patient outcomes in breast cancer.

Methods and Results: Tissue microarrays were used to study p53 IHC expression patterns: expression was defined as extreme positive (EP), extreme negative (EN), and non-extreme (NE; intermediate patterns). Overall survival (OS) was used to define patient outcome. A representative subgroup (n = 30) showing the various p53 immunophenotypes was analysed for TP53 hotspot mutation status (exons 4-9). Extreme expression of any type occurred in 176 of 288 (61%) cases. As compared with NE expression, EP expression was significantly associated (P = 0.039) with poorer OS. In addition, as compared with NE expression, EN expression was associated (P = 0.059) with poorer OS. Combining cases showing either EP or EN expression better predicted OS than either pattern alone (P = 0.028). This combination immunophenotype was significant in univariate but not multivariate analysis. In subgroup analysis, six substitution exon mutations were detected, all corresponding to extreme IHC phenotypes. Five missense mutations corresponded to EP staining, and the nonsense mutation corresponded to EN staining. No mutations were detected in the NE group.

Conclusions: Patients with extreme p53 IHC expression have a worse OS than those with NE expression. Accounting for EN as well as EP expression improves the prognostic impact. Extreme expression positively correlates with nodal stage and histological grade, and negatively with hormone receptor status. Extreme expression may relate to specific mutational status.

The role of secreted frizzled-related protein 2 expression in prostate cancer

AIMS: Improved prostate cancer (PCa)-specific biomarkers are urgently required to distinguish between indolent and aggressive disease, in order to avoid overtreatment. In this study, we investigated the prostatic tissue expression of secreted frizzled-related protein (SFRP)-2.

METHODS AND RESULTS: Following immunohistochemical analysis on PCa tissue microarrays with samples from 216 patients, strong/moderate SFRP-2 expression was observed in epithelial cells of benign prostatic hyperplasia, and negative/weak SFRP-2 expression was observed in the majority of tumour epithelia. However, among Gleason grade 5 carcinomas, 40% showed strong/moderate SFRP-2 expression and 60% showed negative SFRP-2 expression in epithelial cells. Further microscopic evaluation of Gleason grade 5 tumours revealed different morphological patterns, corresponding with differential SFRP-2 expression. The first subgroup (referred to as Type A) appeared to have a morphologically solid growth pattern, whereas the second subgroup (referred to as Type B) appeared to have a more diffuse pattern. Furthermore, 100% (4/4) of Type A patients experienced biochemical recurrence, as compared with 0% (0/6) of Type B patients.

CONCLUSIONS: These results imply: (i) that there is a loss of SFRP-2 expression from benign to malignant prostate glands; and (ii) differential SFRP-2 expression among two possible subgroups of Gleason grade 5 tumours.